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Abstract 

Processive molecular motors which drive the traffic of organelles in cells move in a directed way along cytoskeletal filaments. 
On large time scales, they perform motor walks, i.e., peculiar random walks which arise from the repeated unbinding from 
and rebinding to filaments. Unbound motors perform Brownian motion in the surrounding fluid. In addition, the traffic of 
molecular motors exhibits many cooperative phenomena. In particular, it faces similar problems as the traffic on streets such 
as the occurrence of traffic jams and the coordination of (two-way) traffic. These issues are studied here theoretically using 
lattice models. 
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1. Introduction 

The idea of constructing nanometer-sized devices 
and machines has created a lot of excitement during 
the last years. Despite the progress made, the function- 
ality of artificial nano-devices is, however, still rather 
limited. At the same time, more and more biomolecu- 
lar nano-machines have been identified in the cells of 
living beings where they accomplish a huge variety of 
tasks. Many of these molecular motors are now rather 
well studied and were found to work with an amaz- 
ing degree of precision and efficiency as a result of bil- 
lions of years of evolution, [1,2]. In the following, we 
will focus on one class of molecular motors which has 
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been studied quite extensively during the last decade, 
namely processive cytoskeletal motors which drive the 
traffic of vesicles and organelles within cells. These mo- 
tors hydrolyze adenosinetriphosphate (ATP) and con- 
vert the free energy from this chemical reaction into 
directed movements along filaments of the cytoskele- 
ton. This class of motors contains kinesins and dyneins, 
which move along microtubules, and certain myosins, 
which move along actin filaments. These motors walk 
along the filaments by performing discrete steps with 
a step size which corresponds to the repeat distance of 
the filament, 8 nm for kinesins and 36 nm for myosin 
V. They are called processive if they make many steps 
while staying in contact with the filament. 

From a physical point of view, much of the interest in 
molecular motors is due to the fact that the difference 
in size compared to macroscopic engines implies also 
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conceptual differences. The typical energy of macro- 
scopic motors is much larger than the thermal energy, 
feT, while the typical energies of molecular motors 
are of the order of ksT. For example, the hydrolysis of 
ATP releases about 20 k^T. On the one hand, molec- 
ular motors have to cope with perturbations arising 
from thermal fluctuations; on the other hand, the un- 
avoidable presence of noise suggests that evolution has 
created motors which make use of this noise in order 
to generate work or directed movement. 

From a technological viewpoint, the amazing prop- 
erties of single biological motor molecules and the com- 
plexity of the systems into which they are integrated in 
the cell provide inspiration for the design of artificial 
nanoscale transport systems [3,4]. 

In this article, we discuss several theoretical aspects 
of the motor movements. In section 2, we start by sum- 
marizing some recent experimental results and discuss 
the question whether noise-driven mechanisms are used 
by these motors. In section 3, we discuss another effect 
of noise, namely the detachment of motors from their 
tracks due to thermal fluctuations, which leads to pe- 
culiar random walks. We derive the asymptotic behav- 
ior of these random walks using the statistical proper- 
ties of the returns of motors to the filament. Finally, 
in section 4, we summarize our recent studies of traffic 
problems which arise in systems with many molecu- 
lar motors due to their mutual exclusion from binding 
sites of the filaments. These topics are also addressed 
in a recent longer review article [5] . 



2. Active movements of molecular motors 

Molecular motors can be studied outside cells using 
biomimetic model systems. In these experiments, the 
biological complexity is reduced to a minimal number 
of components, namely motors, filaments, and ATP. 
An example is shown in Fig. 1. By these experiments, 
one can observe movements of single motor molecules 
and measure transport properties such as velocities, 
step sizes, and forces [1,2]. 

On the one hand, these experiments provide insight 
into the motor mechanisms. A major breakthrough was 
to resolve the discrete steps of the motors and to mea- 
sure the step size of the motor which corresponds to 
the repeat distance of the filament. This has first been 




Fig. 1. The 'bead assay' constitutes a biomimetic model 
system: A molecular motor transports a (glass or latex) 
bead along a filament which is immobilized on a surface. 

achieved for kinesin using an optical tweezers setup [6] . 
More recently, it has been shown that kinesin [7,8] and 
also myosin V [9] move in a hand-over-hand way, i.e., 
that the two heads of the dimer step forward in an al- 
ternating fashion such that the rear head always moves 
in front of the leading head, similar to human walking. 
Further progress is expected from combining mechan- 
ical methods and particle tracking with fluorescence 
techniques. 

On the other hand, using these biomimetic motility 
assays, one can measure the transport properties sys- 
tematically varying external control parameters. Here 
the main focus has been on the velocity as a function of 
the ATP concentration and of the force applied with, 
e.g., optical tweezers to oppose the movements, see, 
e.g., Ref. [10]. Other quantities that have been mea- 
sured are the one-dimensional diffusion coefficient of 
motors bound to filaments or the randomness parame- 
ter and the walking distance before unbinding from the 
filament. These measurements have stimulated a large 
amount of theoretical work, see, e.g. Refs. [11,12,13], 
modeling the walks of motors along filaments in order 
either to fit the experimental data or to find out the 
generic properties of these walks. For example, it turns 
out that the motor velocity as a function of the ATP 
concentration is given by a universal relationship which 
should be valid for many types of motors [14,15]. 

Since nanometer-sized molecular motors have to 
live in a noisy environment, it has soon been specu- 
lated whether these motors exploit the noise to gen- 
erate their directed movements, and various variants 
of ratchet models have been proposed as reviewed in 
Refs. [11,12,13]. In the simplest case, the conforma- 
tional changes associated with the chemical cycle of the 
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Fig. 2. Random walk of a molecular motor: The motor 
performs directed movement along a filament and unbinds 
from it after a certain walking distance. The unbound mo- 
tor diffuses in the surrounding fluid until it rebinds to the 
filament and resumes directed motion. 



motor rectify the one- dimensional Brownian motion 
along the filament. While such a simple mechanism 
is not consistent with the measurements for dimeric 
motors such as conventional kinesin or myosin V, it 
can describe movements of processive monomeric mo- 
tors such as the monomeric kinesin KIF1A [16,17,18] 
(for other processive monomeric motors, see the re- 
view [19]). These monomeric motors exhibit biased, 
but strongly diffusive movements along the filaments 
similar to what one obtains in the simplest ratchet 
models. Interestingly, coupling two such ratchets by 
a spring [20,21] leads to a driving mechanism which 
is independent of the diffusion along the filament. 
Similarly, dimerization of monomeric kinesin results 
in a higher velocity and smaller diffusion coefficient 
of the filament-bound motors [17]. These results sug- 
gest that the dimerization of motors has two effects: 
It allows the motor to stay bound to the filament for 
many chemical cycles - in contrast to many (therefore 
unprocessive) monomeric motors - and in addition, it 
allows the motors to move by a more efficient mecha- 
nism which, in contrast to the motility of monomers, 
does not rely on diffusion along the filament. 



3. Random walks arising from many diffusional 
encounters with filaments 

3.1. Motor walks in open compartments 

The fact that nanometer-sized motor molecules work 
in a noisy environment has another consequence: Even 
processive motors do not move along a filament forever, 
but unbind from it after a certain binding time (which 
corresponds to a typical walking distance) , because the 
binding energy is finite and can be overcome by ther- 
mal fluctuations. For a single kinesin or myosin V mo- 
tor, the binding times and walking distances are of the 
order of 1 s and 1 /xm, respectively. Much longer walk- 
ing distances can be obtained if several motors form a 
complex or if a cargo is transported by a larger num- 
ber of motors. Note that the unbinding may also have 
a biological function by allowing the motors to diffuse 
around obstacles on the filament. 

Unbound motors perform simple Brownian motion 
until they rebind to the same or another filament. On 
large time scales, the combination of active directed 
movements along filaments and non-directed Brown- 
ian motion leads to peculiar random walks, called mo- 
tor walks in the following, which consist of alternating 
sequences of the two types of movements as sketched 
in Fig. 2 [22,23]. 

In order to determine the effective transport proper- 
ties of these motor walks, we have studied several sim- 
ple arrangements of filaments embedded in compart- 
ments of various geometries. A particularly simple but 
intriguing case consists of a single filament and a set 
of confining walls, which restrict the diffusion of un- 
bound motors. In the simplest case, there are no con- 
fining walls and the unbound motors can diffuse freely 
in the full three-dimensional space (similar behavior is 
obtained for a half space geometry which is more eas- 
ily accessible to experiments) . By placing the filament 
in a quasi two-dimensional slab or in a cylindrical tube 
(geometries, which are also accessible to in vitro exper- 
iments), diffusion can be restricted along one or two 
dimensions perpendicular to the filament. 

We have studied the motor walks by mapping them 
to random walks on a lattice [23] . A line of lattice sites 
represents the filament. Motors at these sites perform 
a biased random walk and move predominantly into 
one direction, which we choose to be the positive x di- 
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rection. In a discrete-time description, they move with 
a small probability e/2d to each of the adjacent non- 
filament sites and thus unbind from the filament. At 
the non-filament sites the motors perform simple sym- 
metric random walks and move to each neighbor site 
with probability l/2d (d denotes the spatial dimen- 
sion) and rebind to the filament with probability 7r ad 
when they reach again a filament site. Confining walls 
are implemented as repulsive boundaries, at which all 
attempted movements towards the walls are rejected. 

We have used scaling arguments, computer simula- 
tions, and exact solutions of the master equations to 
study the drift and diffusion behavior arising from the 
motor walks [23,24,25,26]. The motor walks exhibit 
anomalous drift behavior and strongly enhanced diffu- 
sion parallel to the filament due to the repeated binding 
and unbinding. At large times, motors move with an ef- 
fective velocity given by fbPb, where Vb is the velocity 
of the bound motor and Pb is the probability that the 
motor is bound to the filament. In the tube geometry, 
Pb is time-independent at large times and given by the 
equilibrium of binding/unbinding and diffusion per- 
pendicular to the filament. In the slab and half-space 
geometries as well as in two- and three-dimensional sys- 
tems without confining walls, a steady state cannot be 
reached at any time, because motors can rebind to the 
filament after arbitrarily large excursions, and the typ- 
ical size of these excursions, that contribute to the aver- 
age behavior, increases with time t. Therefore, Pb and 
the effective velocity are time-dependent in these cases, 
namely Pb(t) ~ t~ d ±/ 2 f or compartments with d± 
dimensions of unconfined diffusion (for d-dimensional 
systems without confining walls, we have d± = d — 1). 
The time-dependent effective velocity implies that the 
average displacement of the motors grows sublinearly. 
In the effectively two-dimensional slab geometry (as 
well as on a two-dimensional lattice without confining 
walls), the displacement behaves as x(t) ~ \/i at large 
times, and in the half space (or full three-dimensional 
space), it is given by x(t) ~ \nt, see Fig. 3. 

In the following, we will give a simplified descrip- 
tion which explains these features and relates them to 
known results from the theory of random walks. 
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Fig. 3. Average displacement function of time 

t for the motor walks as shown in Fig. 2 in two- and 
three-dimensional systems without confining walls. The dis- 
placement grows as x ~ t 1 / 2 and x ~ lnt for d = 2 and 
d = 3, respectively. 

3.2. Asymptotics and return to the filament 

For simplicity, we discuss the case of a single in- 
finitely long filament embedded into d-dimensional 
space with d = 2 or d = 3. Motor particles binding to 
this filament walk along it with velocity t>b until they 
unbind. Unbinding occurs with a rate ~ e, so that the 
motors perform straight movements over the walking 
distance Ax h ~ v h /e. The diffusive excursions be- 
tween two bound walks bring the motors back to the 
filament. (For simplicity, we take the sticking proba- 
bility 7r a( j to be one. If 7r a( j < 1, the motor has to return 
to the filament I/7r ad times before rebinding.) The 
distribution ip(r) of the excursion times r is therefore 
given by the distribution of the return times of a ran- 
dom walker's return to a line in d dimensions or, if we 
consider the projection into the plane perpendicular 
to the filament, of the return to the origin in d — 1 
dimensions. This is a classical problem in the theory 
of random walks which was solved by Polya in 1921 
and has lead to the remarkable result that the return 
of a random walker to the origin is certain on one- and 
two-dimensional lattices, but not in three dimensions 
[27]. For the molecular motors, this implies that they 
will return to the filament with certainty. 

If we are only interested in the movement parallel 
to the filament, we can consider the excursions away 
from the filament as periods of rest. In addition, for 
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large times, the duration of an excursion is typically 
much longer than the time the motor is bound to the 
filament; therefore we can consider the walks along the 
filament as effectively instantaneous steps of size x s ~ 
Axb ~ fb/e '•"ad < 1, the effective step size is x a ~ 
«b7i"ad/e- The motor walks are then described by con- 
tinuous time random walks with the dwell time distri- 
bution V( r ) given by the distribution of the excursion 
times. For this type of random walks, solutions can 
be obtained using Fourier-Laplace transforms [28]. A 
short summary of this method is given in appendix A. 
In that way, one finds that the Laplace transforms of 
the first two moments of the position of this random 
walker are given by 



x{s) = 



and 
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where tf>(s) is the Laplace transform of the waiting time 
distribution ip(s) = / °° dt tp(t)e~ st . In general, {x s } 
and (is) are the moments of the step size distribution; 
for the motor walks, they are given by the active walks 
along the filament and we take them to be given by 
x s = «b/e and x 2 , respectively. From these relations, 
the average displacement and the dispersion of the mo- 
tor walks for large times can be obtained by inverting 
the Laplace transform. In particular, these relations 
imply that the asymptotic displacement of the motors 
is given by the large-time (or small s) behavior of the 
distribution of return times to the origin in d — 1 di- 
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Normal drift behavior with x(t) ~ t is obtained, as 
long as ip(r) has a finite mean value, f , and thus ip(s) w 
I — fs for small s. If, however, ip(r) decays slower than 
~ t~ 2 at large r, the mean dwell time diverges (which 
implies a divergence of [ip(s) — l]/s for small s), and 
anomalous drift is obtained. The latter behavior occurs 
in our case (where the waiting times are given by return 
times to the filament) if the diffusion away from the 
filament is not restricted. 

In the two-dimensional case, the return time distri- 
bution behaves as iP(t) ^ l/(2 v / 7rr 3 / 2 ) for large r or 
tp{s) « 1 — ^/s for small s as shown in appendix B. 
Inserting this into Eq. (1) and inverting the Laplace 
transform, we obtain 
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for small s and large t, respectively. Similarly, in the 
three-dimensional case, the return time distribution 
is ip{r) w 27r/(3rln 2 r) for large r or ip(s) w 1 — 
27r/(31ns~ 1 ) for small s, see again appendix B, which 
leads to 



X(t): 



3x s 



hit 



3«b 

27T6 



Int. 



(5) 



Likewise, we can obtain the dispersion Ax 2 (t) of 
the motors from the second moment of the distribu- 
tion arising from the encounters with filaments. Eq. (2) 
leads to Ax 2 (t) w (v b /e)H + \t and Ax 2 (t) « 

4^r(^bA) 2 lnt + |t in two and three dimensions, re- 
spectively, where we have added the contribution due 
to the diffusion of unbound motors parallel to the fila- 
ment. Note that the broadening of the distribution of 
motors due to the encounters with the filament is char- 
acterized by an anomalously high effective diffusion co- 
efficient of the order of (vb /e) 2 in two dimensions, while 
in three dimensions, the leading term is the unbound 
diffusion, but with a large logarithmic correction, again 
of the order (« b /e) 2 . These results agree with the cor- 
responding asymptotic results from the exact solution 
of the full master equations [24,25]. 



4. Traffic phenomena in many-motor systems 

4.1. Traffic jams and density patterns 

Finally, we consider systems with many interacting 
molecular motors. The simplest type of motor-motor 
interaction is simple exclusion which arises from the 
fact that a motor occupies a certain volume and, in 
particular, excludes other motors from the binding site 
of the filament to which it is bound as observed in dec- 
oration experiments, see, e.g., Ref. [29]. Exclusion is 
most important if motors accumulate in certain regions 
along the filaments, where it leads to the formation 
of molecular traffic jams. These interactions are eas- 
ily incorporated into our lattice model by rejecting all 
movements to occupied lattice sites [23,30,31]. These 
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models then represent new variants of exclusion pro- 
cesses or driven lattice gas models, where the active 
or driven movements are localized to the filaments. In 
driven exclusion processes, the state of the system de- 
pends crucially on the boundary conditions, since the 
boundaries determine the current through the system. 

We have studied tube systems with a single filament 
located along the axis of a cylindrical tube and with 
several types of boundary conditions at the left and 
right tube ends [23,30]. We use the convention that the 
active movement along the filament is biased to the 
right. There are several possibilities to build such tube 
systems artificially such as micropipette glass tubes or 
liquid micro channels, but there are also several tubu- 
lar compartments within cells for which these tube sys- 
tems provide simple descriptions. The most prominent 
example for the latter is the axon of a nerve cell; an- 
other example is provided by the hyphae of fungi. 

The simplest situation is given by periodic boundary 
conditions which can be solved exactly [30] . In this case, 
motor particles reaching the right end of the tube sim- 
ply restart their movements at the left end. This leads 
to constant density profiles for both the bound and 
unbound motors in the stationary state. Diffusive cur- 
rents in the radial direction vanish and binding to the 
filament is locally balanced by unbinding. The current 
of motors through the tube is given by J = Vbpb ( 1 — Pb ) 
and the value of the bound density pb is determined by 
the total number of motor particles in the tube. If the 
number of motors within the tube is increased beyond 
an optimal number (where pb = 1/2 and J = vt>/4), 
the current through the tube decreases due to jamming 
of the motors. 

In a closed tube, motors accumulate in front of the 
right end, and a diffusive current of unbound motors 
to the left balances the current along the filament in 
the stationary state [23,31]. If the number of motors in 
the tube is small, the motors are essentially localized 
at the right tube end. Upon increasing the number of 
motors within the tube, a jammed region at the right 
tube end builds up, separated from a low density re- 
gion to its left by a rather sharp interface, which pro- 
vides probably the simplest example for active pattern 
formation by molecular motors. The crowded domain 
spreads to the left at higher motor concentrations until 
the filament is uniformly covered by motors and rather 
crowded, see Fig. 4. Such density profiles have recently 
been observed for a kinesin-like motor in fungal hyphae 
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Fig. 4. (a) Motors moving in a closed tube system which 
contains one filament. The motors move with velocity Vb 
along the filament to the right and build up a motor traf- 
fic jam at the right end of the system. In the stationary 
state, the bound motor current is balanced by a current of 
unbound motors to the left arising from the diffusive move- 
ment of the unbound motors which is characterized by the 
diffusion coefficient D u b- (b) Profiles of the bound motor 
density as a function of the coordinate x along the filament 
for various total numbers N of motors within the tube. The 
jammed region becomes longer with increasing N . 

[32] . Let us note that these density patterns exist due 
to continuous consumption of ATP. If the ATP concen- 
tration is not kept constant within the closed compart- 
ment, then after burning all ATP, the density pattern 
will finally become homogeneous. 

Tubes with open boundaries which are coupled to 
reservoirs of motors at both ends, so that motors enter 
the tube at the left end and leave it at the right end, 
exhibit boundary-induced phase transitions [30,33]. 
There are three phases which are distinguished by 
the 'bottleneck' which determines the motor current 
through the tube. This 'bottleneck' can be the left 
boundary, the right boundary or the interior of the 
tube. These three cases corresponds to the low-density 
(LD), high-density (HD) and maximal-current (MC) 
phase as shown in Fig. 5 . If changing the motor den- 
sities in the reservoirs at the boundaries leads to a 
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Fig. 5. Phase diagram for the motor traffic through an 
open tube which is coupled to motor reservoirs at the left 
and right end, characterized by bound motor densities pt, ; n 
and pt, cx , respectively [30] . Three phases - a low density 
(LD) phase, a high density (HD) phase, and a maximal 
current (MC) phase - are distinguished by the value of 
the bulk density p{J in the interior of the tube. In the case 
shown here, bound and unbound motor densities satisfy 
radial equilibrium at the boundaries, for other choices of 
the boundary conditions, the transition lines are shifted. 



change in the 'bottleneck' position, a phase transi- 
tion occurs which can be either discontinuous (LD- 
HD) or continuous (LD-MC and HD-MC). These 
types of phases and transitions are known from the 
one-dimensional asymmetric simple exclusion process 
(ASEP) [34,35] which corresponds to the dynamics 
along the filament in our model without the binding 
and unbinding processes. The presence of the unbound 
state of the motors, however, increases the number 
of possible boundary conditions, and since the phase 
transitions are boundary-induced, the location of the 
transition lines within the phase diagram is quite 
sensitive to that choice. At present, it is difficult to 
see which of these boundary conditions will be the 
simplest to implement experimentally. Nevertheless, 
systems of molecular motors are promising candidates 
for the experimental observation of boundary-induced 
phase transitions. 
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Fig. 6. Spontaneous symmetry breaking in systems with 
two motor species moving into opposite directions along the 
same filament: The total motor current J is zero for weak 
interaction q < q c , where the filament is equally populated 
by both motor species, but non-zero for sufficiently strong 
motor— motor interactions with q > q c , where one motor 
species is essentially excluded from the filament [36]. For 
very strong interaction, the current decreases because the 
filament becomes more and more crowded. 

4.2. Phase transitions in two-way traffic 

Another type of phase transitions occurs in systems 
with two-way traffic of motors: While each motor 
moves either towards the plus- or towards the minus- 
end of the corresponding filament, different types of 
motors move into opposite directions. We have stud- 
ied systems with two species of motors moving into 
opposite directions along the same filament [36]. The 
interactions of these motors are described by a single 
interaction parameter q by which the binding and 
unbinding rates are increased or reduced in order to 
enhance binding and reduce unbinding if a motor of 
the same type is present at a neighbor site and to en- 
hance unbinding and reduce binding if a motor with 
opposite directionality is bound at a neighbor site. 
This type of interaction is suggested by decoration 
experiments. For equal concentrations of both motor 
species and for sufficiently strong interaction with 
q > q c (where q c is the critical value of the interaction 
parameter which depends on the overall motor con- 
centration), there is spontaneous symmetry breaking, 
so that one motor species occupies the filament, while 
the other one is largely excluded from it, see Fig. 6. 



7 



If several filament are aligned in parallel and with the 
same orientation, this symmetry breaking leads to the 
spontaneous formation of traffic lanes for motor traffic 
with opposite directionality. Varying the relative con- 
centration of the two motor species for q > q c leads 
to a discontinuous phase transition with hysteresis, 
similar to the transitions induced by changes of the 
external fields in magnetic systems. In contrast to the 
boundary-induced phase transitions discussed in the 
preceding section, these transitions do not depend on 
the boundaries, but are induced by the binding and 
unbinding dynamics. This implies that they are hardly 
affected by the choice of the boundary conditions. 
They depend, however, on the active movements of 
the motors, and are not found in an equilibrium situ- 
ation with motor velocity Vb = 0, as applies, e.g., to 
systems without ATP. 



5. Summary 

In summary, molecular motors exhibit interesting 
movements on several length scales. Here we have ad- 
dressed two of these length scales, namely the walks 
along filaments which consist typically of ~ 100 steps 
with a step size of the order of I Onm and random walks 
which consist of many such walks along filament in- 
terrupted by periods of diffusion after unbinding from 
the filament. In addition, the presence of many motors 
leads to traffic phenomena such as traffic jams and traf- 
fic lanes, similar to the macroscopic traffic on streets. 
However, unbinding of motors from the filaments due 
to thermal fluctuations (which is a consequence of their 
microscopic size) plays an important role and can help 
to circumvent obstacles and to regulate the traffic. 



Appendix A. Continuous time random walks 

In this appendix, we summarize some results for ran- 
dom walks with a dwell time distribution i[>(r) and a 
step distribution T(x 3 ) which are used in section 3. We 
consider random walkers in one dimension which make 
the first step at time t = starting from the origin, 
x — 0. The probability distribution p(x,t) of such a 
random walk fulfills the recursion relation [37,28] 



p(x,i) = 5>„(x) / Vn(t')*(*-Od*'. (A.l) 

n=0 { 

where p„{x) is the probability density that the walker 
is at position x after the n'th step, ip„(t) is the proba- 
bility that the n'th step occurs at time t, and ty(t) = 
f t °° %1)(t)At is the probability that no step occurred un- 
til time t. The initial conditions are po(x) = 8(x) and 
ipo(t) = 5(t). The solution of this recursion can be ob- 
tained using Fourier-Laplace transforms, which leads 
to 

-<*■"= .iii; Sen < A - 2 > 

for the Fourier-Laplace transform of the probability 
distribution p{x,t) [37] with the Laplace transform of 
the waiting time distribution, tp(s) = / °° dt tp(t)e~ st 
and the Fourier transform of the step distribution 
P(q) = J^dx s P(x a )e^ » I + i(x s )q - {x 2 s )q 2 /2. 
The latter expansion is valid for small q, provided 
that the moments {x™} of the step distribution V(x s ) 
with n = 1, 2 are finite. Using this expansion, we can 
derive expressions for the Laplace transforms of the 
moments of our random walk by expanding p(q, s) as 
given in Eq. (A. 2) in powers of q [38] which leads to 
Eq. (I) from which the asymptotic behavior of the 
time-dependent moments can by obtained via the 
Tauberian theorems, see Ref. [28]. 

Appendix B. Return to the origin 

In this appendix we sketch the derivation of the dis- 
tribution of return times to the origin which determines 
the dwell time distribution tp used in section 3.2. We 
consider the probability /(x, t) that, at time t the ran- 
dom walker is at position x for the first time provided 
it started at the origin at time t = 0. /(0, t) then de- 
termines the distribution of return times to the origin. 
These probabilities satisfy the recursion relation 
t 

p(x, t) = /( x , * - t)p(0, t) + 6 t ,oS x ,o (B.f ) 

T = l 

which states that the walker is at site x at time t, if it 
had been there at any time r < t and returned there 
in time t — r. The last term expresses the initial condi- 
tions. Using Fourier-Laplace transforms, one can de- 
rive an expression for the Laplace transform of /(0, t), 
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the probability that the walker returns to the origin at 
time t [27,28], 



f(0,s) 



1 + s (l + s) 2 J( d )(s) 



(B.2) 



where the j( d \s) are integrals over the momentum 
q which depend on the spatial dimensions. The same 
integrals have to be calculated for the exact solution of 
the motor walks [25] . 

In the one- dimensional case we have 



dq 



1 



(1 — cosg) Vs + s 2 
u 

which leads to /(0, s)«l- \fs for small s and 
f(0,t): 



20Ft3/ 2 



(B.3) 



(B.4) 



for large t. In the two-dimensional case the correspond- 
ing integral is given by 



Z7T ZTT 

[S) -(2^J J s + i(2-cos gi 



(2*) 

3^/m 



(2 — cosgi — cos§2) 



K(m), 



(B.5) 



where K(m) is a complete elliptic integral of the first 
kind and m = 4/(2 + 3s) 2 . Since K(m) behaves as 
K(m) m 2 hi f° r m close to one [39], the return 
time distribution is asymptotically given by /(0, s) m 
1 — 27r/(3 In s _1 ) for small s or 



f(0,t) 



2tt 



3iln 2 i 



(B.6) 



for large t. 

The dwell time distribution V>(t) used in section 3.2 
for the case of a filament within a d- dimensional lattice 
is given by i/>(t) = /(0, t = t), where the expression on 
the right hand side has to be taken in the d± = d — 1 
dimensions. 
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